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Abstract At the surface of the Sun acoustic waves appear to be affected by the 
presence of strong magnetic fields in active regions. We explore the possibility that 
the inclined magnetic field in sunspot penumbrae may convert primarily vertically 
propagating acoustic waves into elliptical motion. We use helioseismic holography 
to measure the modulus and phase of the correlation between incoming acoustic 
waves and the local surface motion within two sunspots. These correlations are 
modeled assuming the surface motion is elliptical, and we explore the properties 
of the elliptical motion on the magnetic field inclination. We also demonstrate that 
the phase shift of the outward propagating waves is opposite to the phase shift of the 
inward propagating waves in stronger, more vertical fields, but similar to the inward 
phase shifts in weaker, more inclined fields. 
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1. Introduction 

Helioseismology uses the observed solar surface acoustic wavefield to construct im- 
ages of the subsurface structure of the Sun. Of particular interest has been the 
three-dimensional (3D) modeling of time-distance (Duvall et at, 1993) observations 
of travel-time shifts to deduce the subsurface structure of active regions (Kosovichev, 
Duvall, and Scherrer, 2000; Zhao and Kosovichev, 2003). Assuming the travel-time 
shifts are due to perturbations in the sound speed below the spot, a general con- 
sensus in the models has emerged consistent with sound-speed reductions (relative 
to surrounding quiet Sun) near the surface (< 4 Mm) and enhancements up to 15 
Mm below sunspots (Kosovichev, Duvall, and Scherrer, 2000; Couvidat, Birch, and 
Kosovichev, 2006). Using ring diagram analysis Basu, Antia, and Bogart (2004) also 
find a lower sound speed immediately below the surface and an increase in the sound 
speed below 7 Mm. 
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A comparison between Fourier-Hankel analysis and time-distance results by Braun 
(1997) first prompted caution in the interpretation of acoustic oscillation signals 
within sunspots. The influences of strong surface magnetic fields have not been 
explicitly included in most helioseismic models of active regions. Lindsey and Braun 
(2005a) and Lindsey and Braun (2005b) have shown that helioseismic phase shifts 
observed with helioseismic holography vanish below a depth of about 5 Mm, when a 
surface ( "showerglass" ) phase shift based on photospheric magnetic flux density, is 
removed from the data. 

Other evidence supports the possibility of strong near-surface contributions to the 
helioseismic phase (or travel-time) shifts. These include the possible contamination 
of surface perturbations into the 3D inversions (e.g. Korzennik, 2006; Couvidat and 
Rajaguru, 2007). It has also been shown that the reduction of p-mode amplitudes in 
magnetic regions can cause travel-time shifts (Rajaguru et al., 2006). The suppression 
of sources of wave excitation within sunspots can also produce measurable shifts 
(Hanasoge et al., 2007). 

Schunker et al., (2005) and Schunker, Braun, and Cally (2007) have found that 
phase shifts obtained from seismic holography in sunspot penumbrae vary with the 
line-of-sight angle (from vertical) as projected into the plane containing the magnetic 
field and the vertical direction. A similar effect has also been noted by Zhao and Koso- 
vichev (2006) with time distance measurements. Schunker, Braun, and Cally (2007) 
find that the effect is dependent upon the strength and/or inclination of the magnetic 
field. In the penumbra the magnetic field strength decreases as the magnetic field 
angle from vertical increases, hence the two properties of the magnetic field cannot be 
extricated. The phase variation with line-of-sight viewing angle is demonstrated most 
substantially at frequencies around 5 mHz with a strong, almost vertical magnetic 
field close to the umbra. Schunker, Braun, and Cally (2007) also find that the total 
variation across all lines-of-sight increases with temporal frequency, particularly in 
the stronger fields in the penumbrae. 

Mode conversion of the acoustic waves in the near surface has been explored 
as the physical cause of the observed absorption of acoustic waves by sunspots. A 
fast acoustic wave, propagating towards the surface from the interior, encounters the 
depth at which the Alfven speed is equal to the sound speed (a « c) which is typically 
close to the surface in a sunspot. Under these conditions it is able to transmit to a 
slow acoustic mode and convert to a fast magnetic mode (Cally, 2005). Crouch and 
Cally (2003) explore the mode conversion in two-dimensions with a uniform inclined 
magnetic field relevant to sunspot penumbrae. The inclination of the magnetic field 
is found to have a significant dependence on the likelihood of conversion and fits 
extremely well with the analysis of Braun (1995) (Cally, Crouch, and Braun, 2003). 
Further work (Cally, 2005; Schunker and Cally, 2006) using ray theory has since 
established that it is the angle between the acoustic wave path and the magnetic 
field (the 'attack angle') which is the crucial factor inducing conversion. With a wide 
attack angle at the awe level there is maximum conversion from a fast acoustic to 
a fast magnetic mode. A fine attack angle encourages transmission to a slow acoustic 
mode, which is guided 'up' the magnetic field lines to observational heights in the 
atmosphere. A consequence of mode conversion may be the observational signature 
of elliptical motion in regions of inclined magnetic field. 

The aim of this paper is to model the observations of two sunspots previously an- 
alyzed by Schunker, Braun, and Cally (2007) in order to determine the properties of 
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velocity ellipses consistent with the data. These models are based on a least-squares- 
fit of the phase and modulus information of the local ingression control correlation. 
We explore the properties of these ellipses, observed with waves at different temporal 
frequencies, as functions of the magnetic field inclination angle. We also examine the 
variation with line-of-sight angle of the phase shifts in the outgoing waves, using the 
local egression control correlation, to assess the relation of the phase shift variations 
between incoming and outgoing waves. 

In the following sections we describe the data (Section 2), give an outline of the 
helioseismic holography technique used (Section 3), describe the results (Section 4) 
and discuss our results in the context of mode conversion (Section 5). 

2. Observations 

As this is a continuation of studies by Schunker et a/., (2005) and Schunker, Braun, 
and Cally (2007) we use the same data, however, we provide a short description here 
to maintain coherence. The Michelson Doppler Imager (MDI) aboard the Solar and 
Heliospheric Observatory (SOHO) (Scherrer et a/., 1995) provides the solar surface 
Doppler velocity information. The Dopplergrams are full disk, have a 60 second 
cadence and a resolution of ss 1.4 Mm per pixel. These full disk Dopplergrams are 
Postel projected and a 512 x 512 pixel extract is taken centered on the active region. 
We analyze two sunspots: the first in AR9026 observed over 10 days from 3rd - 
12th June 2000 with a Carrington longitude (L0) of 75° and latitude (B0) of 20° 
and penumbral boundaries defined by inner and outer radii of 7 Mm and 16 Mm 
respectively. The second is a sunspot in AR9057 observed over 9 days from 24th 
June - 2nd July 2000 with L0 = 158° and B0 = 13°, inner and outer penumbral 
boundaries given by 6 Mm and 13 Mm. The penumbral boundaries are defined to be 
between 50% and 85% of the nearby quiet-Sun continuum intensity. The sunspots in 
these active regions were chosen based on the existence of continual MDI observation 
as they traversed the solar disk, and for their relatively simple magnetic structure 
and evolution. 

Magnetograms from the Imaging Vector Magnetograph (IVM) at the University 
of Hawaii Mees Solar Observatory (Mickey et a/., 1996) provide the orientation and 
strength of the surface magnetic field in the sunspots chosen in AR9026 and AR9057. 
The IVM observations are made over a 28 minute interval: for AR9026 starting at 
18:29 Universal Time (UT) on 5th June 2000 and for AR9057 starting at 16:19 
UT on 28th June 2000. The IVM data reveals an azimuthally spreading magnetic 
field configuration for both sunspots although they are not entirely symmetric. It 
is assumed, supported largely by available line-of-sight magnetograms, that there 
is no significant evolution of the magnetic field in the sunspots during the time of 
observation. Therefore, using only one vector magnetogram for the duration of the 
observation is reasonable. Rotation and scaling are applied to align the IVM data to 
the line-of-sight MDI magnetograms. 

Figure 1 shows a strong correlation between the magnetic field inclination (from 
vertical, 7) and field strength (the strong field is almost vertical whereas highly 
inclined field is relatively weak). In this paper we use the inclination 7 as the pri- 
mary variable, dividing the penumbrae into three regions defined by the values of 7 
(Fig 1). It is understood in this analysis that the magnetic field strength is implicitly 
correlated with the inclination through Figure 1, and that independent dependencies 
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of observables with field strength and inclination are not extracted. The dependence 
of the phase shifts on the line-of-sight viewing angle of the magnetic field is facilitated 
by knowing the full vector magnetic field. 
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Figure 1. Magnetic field strength, \B\, plotted against inclination from vertical, 7, determined 
from IVM vector magnetograms. AR9026 (5th June 2000) (bottom) and AR9057 (28th June 
2000) (top). The different symbols divide the penumbra into roughly equal regions of inclina- 
tion: 7 < 42° (asterisk); 42° < 7 < 66° (diamonds); 7 > 66° (triangles). This corresponds to 
average magnetic field strengths of 1700 G, 1000 G and 600 G for AR9057 and for AR9026 
1900 G, 1400 G and 600 G. In general, the progression from the upper left portion of the 
distribution to the lower right portion represents increasing distance from the centre of the 
spot. 



3. Helioseismic Holography 

Helioseismic holography (Lindsey and Braun, 2000; Braun and Lindsey, 2000) is the 
phase coherent imaging of the solar subsurface based on photospheric acoustic oscil- 
lations. The ingression is an assessment of the observed wavefield, ^(r , t), converging 
to a selected focal point, (r, z, t), and the egression is the time reverse - an assessment 
of waves diverging from that point. In this case we calculate the quantities at the 
surface, z = 0. In practice, the observed wavefield used for the calculation is usually 
an annulus surrounding the chosen focal point. The pupil used here is identical to 
that described by Schunker et al, (2005) and is constructed for the calculations 
with inner radius a = 20.7 Mm and outer radius b = 43.5 Mm, designed to be large 
enough that when the focal point is within the penumbra the area covered by the 
annulus does not include large areas of strong magnetic field. At a frequency of 5 
mHz this selects p-modes with spherical harmonic degree and radial degree between 
I « 450 and I rj 700. The ingression at the surface is given by, 

H-(v,0,t)= [ dVG-(|r-r'|,0,£-O<Kr',t), (1) 

J a<\r-r'\<b 
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G- is the ingression Green's function (Lindsey and Braun, 2000) . The egression, H+ 
is simply the time reverse of 1, i.e. t — t. 

The ingression (1) is correlated with the observed wavefield in the space-frequency 
domain, 

C-{v,v) = {H-{v,v)^{v,v)) Av = \C-\e~' 5<t, ~ , (2) 

where H-(r,i/) and ip(r, v) are the temporal Fourier transforms of H- and ip re- 
spectively, and we have dropped the dependence on depth z. The correlation has a 
modulus |C-| and a phase, 5<j)-- In the frequency spectrum, with the pupil covering 
mostly quiet Sun, the local ingression control correlation simply characterizes how 
the local magnetic photosphere responds acoustically to upcoming waves, prescribed 
by H- , originating within the pupil. 

In Schunker et at, (2005) and Schunker, Braun, and Cally (2007) the phase shift 
caused by the surface perturbations to the incident wave is calculated at various 
frequencies. Here we extend that research and monitor the variation of the correlation 
amplitude within the sunspot penumbral region, which is then combined with the 
phase information to form estimates of the surface velocity ellipse. We then go on to 
examine the phase of the local egression control correlation, 5<f>+ = Arg[C+(r, v)], 
with the line-of-sight angle in the penumbra of sunspots in AR9057 and AR9026. 



4. Elliptical Representation of Surface Velocities 

For better statistics multiple days of observation of each sunspot are combined and 
a least-squares-fit of the observations allows an estimation of the velocity ellipse. To 
calculate the velocity ellipse the ingression correlation (eqn. 2) is used as a proxy 
for the local velocity. The modulus and phase of the correlation vary with respect 
to the line-of-sight angle and these variations are modeled to construct an ellipse 
representing the surface velocity associated with a particular magnetic field element. 

A smeared ingression 'flat field' takes out undesired contributions to the corre- 
lation modulus due to the temporal or spatial variations of the ingression, due to 
the presence of magnetic regions in the pupil. This is achieved by dividing |C_| by 
the Gaussian smear of the root-mean-square of the ingression to get |C_|^; a4 . The 
resulting correlations are normalized so that the ingression correlation modulus in 
the nearby quiet Sun has a value of cos(£), where £ is the heliocentric angle of the 
quiet region from disk center. A dependence of the quiet Sun (root-mean-squared) 
amplitude with cos(£) is expected for predominantly vertically oscillating p-modes. 
The normalization is achieved by dividing \C-\fi a t by the quiet Sun average divided 
by the cosine of the heliocentric angle, i.e. 

| C- 1 norm = \C-\fl a t 7T" • (3) 

{\G-\)qs 

The normalization corrects for variations in the modulus due to duty cycle variations, 
foreshortening, and other effects which cause undesired variations in the modulus 
from day to day. These procedures assume that these factors have the same relative 
effect on the (desired) correlations in the penumbra as they do to the quiet Sun. It 
is difficult to assess the validity of this, and so it is used as a reasonable working 
assumption subject to some caution. 
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We use the same angle O p as defined in Schunker et al, (2005), which is the 
angle between the projection of the line-of-sight vector onto the plane containing 
the magnetic field vector and the radial vector. It is a measure of the angle that the 
magnetic field is being viewed from. We now simply refer to |C_| raorTO as \C- \ and 
see how it varies with 9 V . 
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Figure 2. The modulus of the correlation, \C-\ (left column), and the phase, &4>—, in the 
penumbra of AR9026 at 3 mHz for all days of observation are plotted against projected angle 9 P 
for different values of magnetic field inclinations as indicated. The top panel (a) shows 7 < 42°, 
where the mean field strength is (B) = 1900 G, the middle panel (b) shows 42° < 7 < 66°, 
where <B) = 1400 G, and the bottom panel (c) shows 7 > 66°, where <B> = 600 G. The 
horizontal dashed lines indicate the mean value of |C— | for each panel. The error bars indicate 



a fit for all the displayed data; the dotted line is a fit for the data from 3rd - 7th June 2000; 
the dashed line is a fit for data from 8th - 12th June 2000. 



Figures 2 to 4 and Figures 6 to 8 show the variation of \C- | with P for AR9026 and 
AR9057 in the left columns at 3, 4 and 5mHz in the same three bins of inclination 
shown in Figure 1. The right column shows 5(j>—. The variations of 5<j)— with 9 P 
have been the subject of our previous analyses (Schunker et al., , 2005; Schunker, 
Braun, and Cally, 2007). In the absence of magnetic effects, we expect the local 
oscillatory wave field, as assessed by the ingression correlation, to be consistent with 
purely vertical motion. This would predict a dependence of |C_| on cos(0 p ) and no 
dependence of 5(f>— on cos(6 p ). Departures from these expectation are clearly visible 



thesispaper_db.tex; 22/04/2008; 10:50; p. 6 



Surface Velocities in Magnetic Fields 



?:8 



5 :8 



(0)7 < 42* 




• 

(b)42' < 7 < 66' 


• 

i 


• 

(c)7 > 66' 


— i 1 — i — i 1 — 



-75 -25 25 75 -75 -25 25 75 
B, (deg) 6, (deg) 

Figure 3. Same as Figure 2 except at 4 mHz. 



in Figs 2-4 and 6 - 8. A rough understanding of these results can be gained by 
noting that the net variation in the phase shift with 6 P is related to the eccentricity 
of the ellipse, while the value of 9 P for maximum \C- \ determines the orientation of 
the semi-major axis. We perform a least-squares- fit to the observed \C-\ and &4>- 
to determine the elliptical motion consistent with the data. The best-fit ellipses are 
shown in Figures 5 and 9. The moduli and phase of C— as determined from the 
fits are plotted (as solid lines) with the data points in Figs 2-4 and 6-8. Separate 
fits were performed for independent 5-day subsets of the data (also shown in the 
Figures by the dotted and dashed lines). There are no obvious systematic differences 
in the fits over time. The orientations and eccentricities are highly consistent across 
all frequency bands in both sunspots (Figure 5 and Figure 9). Some systematic 
differences between the two spots are evident. For AR9026, the ellipses are aligned 
slightly towards the magnetic field direction at high field inclination but 'swing' over 
as the inclination becomes smaller (and the magnetic field stronger). For AR9057, 
the motion is nearly vertical at low field inclination, but also tilts away from the field 
in the stronger, more vertical, fields. In both spots, the eccentricity of the ellipses 
increases with decreasing field strength or increasing inclination. 
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4.1. ELLIPSE PARAMETERS 

We define a deviation angle as the angle between the magnetic field vector and surface 
velocity ellipse semi-major axis (5 = 7-/3, where f3 is the inclination of the semi- 
major axis from vertical). We show the variation of the deviation angle, length of 
the semi-major axis, and eccentricity with magnetic field strength and/or inclination 
in Figure 10. Also shown is the phase difference between the fits at 9 P = —60° and 
Op = +60° . This quantity is generally inversely correlated to the ellipse eccentricity, 
but is a more direct measure of the total variation of observed phase shift for a 
specific penumbral region (Schunker et al, 2005). In Figure 10, data from AR9026 is 
represented by an asterisk, AR9057 by a diamond and the frequencies are color-coded 
as following :- 5 mHz is black, 4 mHz is purple and 3 mHz is red. 

In general, the trends shown in Figure 10, namely an increase in the deviation 
angle, eccentricity, and semi-major axis length, and a decrease in the phase variation, 
with increasing inclination are observed in both sunspots and at all frequencies. There 
are some deviations from this. For example, at low inclinations, it is observed that the 
eccenctricity (phase variation) decreases (increases) with frequency. Note that the 
semi-major axis length is an indication of the total wave amplitude in the magnetic 
region. The trend observed in the lower-left panel of Figure 10 is consistent with a 
reduction in wave amplitude related to the field strength. 
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Figure 5. Lcast-squares-fit surface velocity ellipses as given by the phase and amplitude of the 
local ingrcssion control correlation in the penumbra of AR9026. The top row is when 7 > 66° , 
the middle row when 42° < 7 < 66° and the bottom row is when 7 < 42°. The 7 listed in the 
plot is the angle that the magnetic field vector is drawn at. (3 is the inclination angle of the 
semi-major axis of the velocity ellipse. The left column is at frequencies of 3 mHz, the middle 
column at 4 mHz and the right column at 5 mHz. 

5. Local Egression Control Correlation 

In this section we examine the phase of the local egression control correlation, 

C+{v,v) = {H + (t,v)-4>*{t,v)) Av = IC+le- 1 4 *+. (4) 

Since the egression is simply the time reverse of the ingression, we might expect to 
see a reversal of the phase change compared to the ingression phases. 

Using all the days' of data, the egression correlation phase is plotted against 9 P in 
Figures 11 and 12, along with the fits for the phase variation due to elliptical motion. 
We see similar trends in both sunspots, AR9026 and AR9787. Figure 13 shows the 
phase difference between the least-squares-fit correlation phase at 9 P = —60° and 
at P = +60° of the ingression plotted against that of the egression. The colors 
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Figure 6. The modulus of the correlation, \C— | (left column) and phase, Sip— (right column), 
in the penumbra of AR9057 at 3 mHz for all days of observation are plotted against projected 
angle 9 P for different values of magnetic field inclinations as indicated. The top panel (a) 
shows 7 < 42°, where the mean field strength is (B) = 1700 G, the middle panel (b) shows 
42° < 7 < 66°, where (B) = 1000 G, and the bottom panel (c) shows 7 > 66°, where 
(B) = 600 G. The horizontal dashed lines indicate the mean value of \C— \ for each panel. The 
error bars indicate the standard deviation of the mean over bins of 20 measurements in 8 P . 
The solid line of fit is a fit for all the displayed data; the dotted line is a fit for the data from 
24th - 28th June 2000; the dashed line is a fit for data from 29th June - 2nd July 2000. 



and symbols represent the same as before: AR9026 is represented by an asterisk, 
AR9057 by a diamond and the frequencies are 5 mHz (black), 4 mHz (purple) and 
3 mHz (red) in addition the size of the symbols represent the average inclination 
from vertical. The solid line has a slope of —1. There is a reverse behaviour of the 
ingression, compared to the previous egression results, present for all frequencies at 
most magnetic field inclinations. The exception is when the field is highly inclined, 
where we observe a trend in the same sense as the ingression. This is unexpected 
and warrants further study. 



6. Discussion 



Mode Conversion predicts (among other things (Cally, 2007)) that when the attack 
angle is small most of the observable energy will be in the slow acoustic mode, 
and when the attack angle is large most of the observable energy will be in the fast 
magnetic mode. In this case we will be seeing the line-of-sight effect of a combination 
of waves coming from all directions impinging on magnetic field with a particular 
orientation. The observations by MDI consist of line-of-sight Doppler signatures of 
the surface motion which are presumably caused mainly by pressure perturbations. 
This means that we would expect to be observing only the slow acoustic mode, 
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Figure 7. Same as for Figure 6 except at 4 mHz. 



however it is possible that we are observing a combination of the acoustic and 
magnetic modes. 

The dependence of the phase shift of the observed ingression correlations with 
azimuthal angle around sunspot penumbrae, as viewed from different observational 
vantages, shows that the incoming phase shifts must be (at least partly) photospheric 
in origin and are influenced by the presence of inclined magnetic fields (Schunker 
et at, 2005). Analysis of the variation of both the amplitude and phase of the 
surface velocities provides an opportunity to characterize the magnetically influenced 
acoustic signature as an ellipse with properties determined by the magnetic fields. We 
found that the ellipses are either nearly vertical (for weaker, more inclined fields) or 
generally directed away from the magnetic field direction (for stronger, more vertical 
fields). Largely consistent results for two active regions, AR9026 and AR9057, are 
found. Some properties of the surface ellipses, e.g. their inclinations, are different 
for the two sunspots. Some of this variation may be due to differences in the field 
properties. For example, the field in AR9057 is on average ~ 15% weaker than 
AR9026. 

Fits of the elliptical motion in Figures 5 and 9 depend critically on the correlation 
modulus which is prone to systematic uncertainties. But the trend is that a stronger, 
less inclined magnetic field produces elliptical motion with smaller amplitude, eccen- 
tricity, and deviation angle, and a larger inclination from vertical. These trends exist 
for both spots and, largely, at all frequencies. The shorter semi-major axis at strong 
magnetic field strengths is consistent with previous knowledge of surface acoustic 
amplitude suppression in magnetic fields, ft is curious to note, however, that at 4 
mHz the amplitude is consistently smaller than even 3 mHz. 

These are the first results to estimate the behaviour of the surface velocity ellipse 
at the photosphere within sunspots. The results do not immediately suggest an 
observation of the slow wave as shown by Cally (2005) or Schunker and Cally (2006), 
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but are consistent with their expectations of the behaviour of slow waves at this 
height in the atmosphere. Mode conversion theory states the alignment is dependent 
upon a 2 /c 2 which at the observational heights of ~ 200km of the atmosphere may not 
be large enough to invoke alignment. Since the ray analysis is somewhat unrealistic 
we would expect to observe a combination of fast and slow waves at the surface, which 
will contribute to a clouded view of the surface velocities. Rajaguru et al, (2007) are 
currently exploring the possibility that these apparent surface effects are due to the 
changes in the radiative transfer within active regions and the formation height of the 
observational Ni 678 nm line. This explanation requires an absorption mechanism, or 
else some other means of producing a difference between the amplitudes of upward 
and downward propagating waves. Thus mode conversion may still be important 
in this proposed mechanism. A test of the mechanism proposed by Rajaguru et 
al. (2007) would be to repeat the observations performed here in a magnetically 
insensitive line, where the proposed radiative transfer effects would not be present. In 
terms of mode conversion, it is suggested that the main effect occurs along the bright 
radial filaments of the interlocking comb structure as presented in the penumbral 
models of Weiss et al, (2004). However, observational helioseismic spatial resolution 
cannot currently resolve this. 

This is also the first time that the variation of the phase of the local egression 
correlation has been analysed in the penumbra. It is curious that the egression 
correlation shows a reverse dependence when the magnetic field is weak and highly 
inclined. This is evidence of a reverse ingression dependence on the line-of-sight, but 
further investigation is required to understand the behavior at high frequencies in 
the weaker, more inclined fields. 
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Figure 9. Lcast-squares-fit surface velocity ellipses as given by the phase and amplitude of the 
local ingrcssion control correlation in the penumbra of AR9057. The top row is when 7 > 66° , 
the middle row when 42° < 7 < 66° and the bottom row is when 7 < 42°. The 7 listed in the 
plot is the angle that the magnetic field vector is drawn at. (3 is the inclination angle of the 
semi-major axis of the velocity ellipse. The left column is at frequencies of 3 mHz, the middle 
column at 4 mHz and the right column at 5 mHz. 
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Figure 11. The 3 mHz (left column), 4 mHz (middle column) and 5 mHz (right column) 
egression correlation phase (<5</>+) vs. 8 P within the penumbra of sunspot AR9026 for different 
values of magnetic field inclination as indicated. The three rows represent different portions of 
the penumbra as shown in Figure 1(a). The top panel row (a) shows 7 < 42°. where the mean 
field strength is (B) = 1900 G, the middle row (b) shows 42° < 7 < 66°, where (B) = 1400 
G, and the bottom row (c) shows 7 > 66°, where (B) = 600 G. The horizontal dashed lines 
indicate the mean value of 5<f+ f° r each panel. 
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Figure 12. The 3 mHz (left column), 4 mHz (middle column) and 5 mHz (left column) 
egression correlation phase (<5</>+) vs. 8 P within the penumbra of sunspot AR9057 for different 
values of magnetic field strength as indicated. The three different panels represent different 
portions of the penumbra, similar to Figure 1(b). The top panel (a) shows 7 < 42° where 
the mean field strength is (B) = 1700 G the middle panel (b) shows 42° < 7 < 66°, and 
(B) = 1000 G and the bottom panel (c) shows 7 > 66° and (B) = 600 G. The horizontal 
dashed lines indicate the mean value of <5<yi_|_ for each panel. 
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Figure 13. The phase difference between the least-squares-fit correlation phase at 8 P = —60° 
and at 8 P = +60° of the ingrcssion plotted against that of the egression. The solid line is a 
line of slope —1. AR9026 is represented by an asterisk, AR9057 by a diamond and the colors 
indicate the following frequencies: 5 mHz (black), 4 mHz (purple) and 3 mHz (red). The largest 
symbols represent 7 > 66° and the smallest symbols 7 < 42° . 
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